Thermal phase diagrams of columnar liquid crystals 
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In order to understand the possible sequence of transitions from the disordered columnar phase to the 
helical phase in hexa(hexylthio)triphenylene (HHTT), we study a three-dimensional planar model 
with octupolar interactions inscribed on a triangular lattice of columns. We obtain thermal phase 
diagrams using a mean-field approximation and Monte Carlo simulations. These two approaches 
give similar results, namely, in the quasi one-dimensional regime, as the temperature is lowered, the 
columns order with a linear polarization, whereas helical phases develop at lower temperatures. The 
helicity patterns of the helical phases are determined by the exact nature of the frustration in the 
system, itself related to the octupolar nature of the molecules. 
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I. INTRODUCTION 

The study of phase transitions in columnar liquid crys- 
tals presents a fundamental interest in that these 
materials combine, aside from the vast phenomenology 
of soft matter, many features at the origin of impor- 
tant phenomena of the solid state, namely a relatively 
strong elastic anisotropy in the direction of the columns, 
a geometrical frustration of the intermolccular interac- 
tion coming from the triangular nature of the lattice 
of columns, and finally discoid molecules with nontriv- 
ial point-group symmetry. The present study is based 
on hexa(hexylthio)triphenylcnc (HHTT), a compound 
formed of a rigid core of aromatic cycles, giving a dis- 
coid shape to the molecule, and of six flexible hydrocar- 
bon chains, responsible for its characteristic thermotropic 
character. HHTT is the only molecule from the tripheny- 
lene derivatives to show two distinct columnar phases. 
Indeed, as the temperature decreases, the sequence of 
phases is the following: /, an isotropic liquid; D^d, a dis- 
ordered columnar phase; H, a helically ordered columnar 
phase and finally K, a monoclinic crystal. These phases 
were identified by X-rays measurements on powders H 
and freely suspended strands of HHTT. 

These X-rays results are best interpreted by conclud- 
ing that the D hd phase of HHTT (70°C< T < 93°C) 
has long-range positional order in the plane perpendicu- 
lar to the columns: columns are located on a triangular 
lattice. Within a column, short range (liquid) positional 
order is realized. The columns slide freely one against the 
other. The H phase (62°C< T < 70°C) has in-column 
long-range positional and orientational helical order, or 
quasi long-range order, as proposed in 0. In this last 
phase, two neighboring molecules in a single column are 
separated on the average by a distance d\\ — 3.6A and 
rotated from each other by an angle a « 45°, constant 
on the whole temperature interval of the phase. In or- 
der to reduce the frustration associated with the trian- 
gular geometry of the lattice, the lattice reorganizes it- 
self in a superlattice \/3 x \/3i?30 o : one third of the 
columns have a vertical offset of half the inter- molecular 
distance (dii/2). The displaced columns have an opposite 



helicity: a w —45°, instead of +45° for the undisplaced 
columns. If it were not for the very high value of the mean 
square displacement of the molecules in the direction of 
the columns, the H phase would seem very similar to a 
crystalline phase. In that sense, the exact mechanism of 
the Dhd <-> H transition is still an open question. In an 
effort towards elucitating the mechanism, we study the 
orientational ordering assuming from the start that the 
HHTT molecules already occupy well defined positions. 

We study the thermal phase diagram of a model Hamil- 
tonian by means of a Landau free energy functional in 
a mean-field approximation and of Monte Carlo simu- 
lations on finite size lattice. Previous work has been 
done on the ground state of a related model J^] and on 
thermal phase diagrams for a two-dimensional model of 
uniform columns p]lCf|. Our analysis confirms that, as 
previously seen at T = 0, the octupolar G coupling O] 
is determinant for obtaining the helicity configuration of 
the columns at any temperature. It also shows that a 
diversity of phases survives at T ^ 0. For weak trans- 
verse couplings, the model produces the expected low- 
temperature helical phases, but also suggests that some 
linearly polarized phases could exist at higher tempera- 
tures. 

In the following section, the model Hamiltonian is pre- 
sented with an emphasis on the inter-columnar inter- 
molecular interactions. In Sec. Ill, using mean-field ap- 
proaches, the locus of points for the second order phase 
transitions from the disordered phase to an ordred phase 
is obtained. It is followed by a determination of the ther- 
mal phase diagrams and characterisation of the differ- 
ent phases in terms of the helicity pattern and relative 
phases. In Sec. IV, the thermal phase diagrams are ob- 
tained using Monte Carlo simulations on finite size sys- 
tems in conjunction with the spiraling algorithm. Finally, 
Sec. V discusses the results and arrives at general con- 
clusions. 
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II. MODEL HAMILTONIAN 

As indicated it the introduction, the main purpose of 
our calculation is to elucidate the role played by the an- 
gular degrees of freedom. Even though this may be ques- 
tionable formally, we will freeze the positional degrees 
of freedom to render the problem tractable. Despite the 
fact that the exact positional order in the columns is 
not known, we then assume that the molecules lie on a 
three-dimensional triangular lattice of ordered columns 
(with one third of the columns being displaced). This 
simplification is certainly valid in the H phase, even if 
only quasi long-range positional order existed since or- 
der would then be maintained over many inter-molecular 
distances. Accordingly, we will use site indices to 

identify the unit cell of the three columns of the two- 
dimensional V3 x V3-R30 superlattice of columns and [i 
to identify the column: \i = 1 and 2 label the undisplaced 
columns and fi = 3 labels the columns offset by d\\/2. Fi- 
nally, an index k identifies the sites a molecule occupies 
in the (p,i,j) column. We formally write (i,j,k) as m, 
an index labelling a plaquette of three molecules. 
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FIG. 1. a) Two-dimensional triangular lattice of columns. 
Empty dots represent displaced columns (/i = 3) and 
filled dots, undisplaced columns (/x = 1 and 2). Dashed 
lines represent the elementary cell of the super-lattice, b) 
Three-dimensional illustration of the couplings for three 
columns. Dotted lines represent Active sites and interactions 
for an undisplaced column 3. 

Within the framework of this plastic state model, each 
molecule (/i, m) has a well-defined position r Mm and an 
orientation labelled 0^ m , admitting from the beginning 
that its plane is perpendicular to the direction of the 
columns. An orientational disorder may mimick an ef- 
fective Dhd phase, which however would possess orien- 
tational and positional disorders along the columns. We 
did not explicitly consider the shape and flexibility of the 
tails, which would vary with the temperature. Recall, 
however, that it has been suggested that the stiff- 
ening of the tails may be responsible for the Dhd —> H 
transition. The above considerations are implicitly inte- 
grated out as weakly temperature-dependent renormal- 
ization effects of the inter-columnar interactions, allowing 



us to use, to a good approximation, an effective model 
with athermal values of the interaction parameters. The 
resulting thermal phase diagrams will represent a some- 
what distorted version of the true temperature depen- 
dence. 

The last thing to consider in the model is that the 
molecules are not exactly disk-like: they have a chi- 
rality associated with the alternate arrangement of the 
aliphatic tails. Indeed, con forma tional analysis on corn- 
shows that this 



pounds similar to HHTT |12 
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pro- 



peller blade" configuration is the ground state of a single 
molecule and of two stacked molecules, one on top of the 
other. 



A. Intra-columnar interactions 

In their ground state, two stacked HHTT molecules 
minimize their conformation energy by allowing an an- 
gular shift a between the two molecules [fiJlO!. An in- 
trinsic chiral model represented by the Hamiltonian 



Jcos3(6» fc+ i -9 k -a), 



(1) 



would be appropriate to represent this situation since the 
HHTT molecules have D% symmetry. However there is no 
a priori selection between the right-handed or left-handed 
chirality. In order to allow for the two possible chirali- 
ties of each column (the sign of a), we preferably use a 
next-nearest-neighbor model (see p^ , p"5[ ), with compet- 
ing interactions of the form: 

- Ji cos 3(0 fc+ i - 6 k ) - J 2 cos 3(6> fc+2 -0k), (2) 

with J\ > and J2 < 0. In the ground state, for 
4|J2| > \Ji\, the molecules adopt a helical configuration 
with pitch a, given by cos 3a = — J1/4J2. Otherwise, 
the intrinsinc helicity of the columns is zero. The J1/J2 
ratio determines the magnitude of a, but allows opposite 
helicities for different columns or even helicity reversals 
within a column, separating helicity domains. 

For an isolated molecule of HHTT at finite tempera- 
ture, it is unclear that chirality is a well defined property: 
conformational analysis jl2| shows that the energy bar- 
rier between opposite chirality configurations is compara- 
ble to the thermal energy in the H and Dhd phases. Thus, 
instead of using an additional Ising variable on each site 
to represent the chirality of the molecule, we represent 
the chirality as the result of a collective behavior in the 
effective model (0), submitted to inter-columnar inter- 
actions: every molecule in a particular helicity domain 
has the same chirality, related to the sign of the helicity. 
The inter-columnar interactions are essential to stabilize 
helicity domains at finite temperature. 



B. Inter-columnar interactions 

Given the approximation that each molecule is fully 
described by its orientation 9^ m , we may write its mass 
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density as a multipole expansion 0,^6) . Because of the 
Z?3 point symmetry of the molecule, the first nonzero mo- 
ment is the octupolar moment, which may be represented 
by Qkim, a rank-three tensor (k,l,m = x,y). The only 
interactions that are bilinear in Q as well as invariant 
with respect to the symmetries of the hexagonal lattice 
have a cos3(6> - 6') or cos3(6> + 9') form Q. The inter- 
columnar interaction is then approximated to be 



£ 
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J l3 (q)(S 3 (q)), 



(J) 



where Jij (q) is a 6 x 6 block diagonal matrix constructed 
from the J c and J s couplings of Eq. (0) : 



J(q) = 



J c (q) 






J s (q) 



(8) 



Jcos3(6*p 



Gcos3(# p 



(3) 



The first term is invariant under continuous rotations 
and would be present even if the molecules had lower 
symmetry multipole moments. However, the second term 
is specific to the octupolar character of the molecules and 
has only discrete rotational symmetries. 

To extract the D$ symmetry of the molecules, we 
replace the real orientations 9^ m by angular variables 
4>iim — 3# Mm . The complete hamiltonian of the system 
then reads: 



H 



fiv mn 



Jfim,un contains the intra-columnar interactions: 



(4) 



each 



site is coupled to its first and second intra-column near- 
est neighbors by J\ and Ji as in Eq. (§). The inter- 
columnar couplings are embedded both in J U m,un and 
Each n — 1 site interacts with three [i = 2 



G 



and six /i = 3 neighbors, three upwards and three down- 
wards. The inter-columnar couplings have different val- 
ues: J and G for in-plane molecules (1-2 bonds) and J' 
and G' for out-of-plane molecules (1-3 and 2-3 bonds) 
(see Fig. |l| b). Nevertheless, it is physically justified |J 
to suppose that J' » J and G 1 « G, and for simplicity, 
we assume J' = J and G' = G. J\ is positive and taken 
to be unity (it sets the energy scale). 
With the notations 



c M m = cos(</y m ) and s Mm = sin 
we may re-write the hamiltonian (m as 



0, 



H = 



au mn 



(5) 



(6) 



Because there are three columns in the unit cell, the 
mean-field transverse components of q are zero, as we 
have verified through a detailed calculation. From now 
on, without any ambiguity, we replace q by q, its z com- 
ponent. We also take the intra-column distance between 
two molecules d\\ =1. The J c and J s matrices are then 



J cs (q) 



' J\\(q) 

jC,S 

• J 3i S (?) 



jC,S 

J A q) 

^23 S (<?) 



J 31 s (l)' 

J\\(q) . 



(9) 



with Jii (q) = cos q + J 2 cos 2q, Jf 2 = |(J + G), J| 3 (g) = 
JiM = 3(J + G) cos \q, j; a = |(J - G) and Jf 3 (g) = 
J§i{q) — 3(J — G) cos iq. The displacement of the (X = 3 
columns changes the coordination number from 6 to 3 
and adds a cos \q factor. 



A. Second-order phase transition temperature T c 

The order-disorder continuous phase transitions are re- 
lated to the divergence of the "paramagnetic" suscepti- 
bility x, which in turn is related to the single-site suscep- 
tibility xo = 1/2T, with k B = 1, by the standard RPA 
relation: 



x(q) = xo[a - xoJ(q)Y 



(10) 



The 6x6 matrix between brackets is non-invertible when 
at least one of its eigenvalues is zero. As the temperature 
is lowered, the transition occurs for some q c maximizing 
one of the six eigenvalues of J(q). The corresponding 
eigenvector identifies the configuration involved in the 
transition. The eigenvalue itself is twice the critical tem- 
perature T c . 

J(q) is block-diagonal, and the 6th-order caracteristic 
equation reduces to the two cubic equations given by 



where J n m,vn Jnm,V7i~\~G^ 1 fy ljj vn and J am un — Jum,i>n 

G am ,vn- The reader should note that the G <-* — G trans- 
formation interchanges the c and s variables and is equiv- 



alent to a rotation 



| of the molecules. 



det(J c '*(<z)-j c ' s l)=Q. 
The eigenvalues are 



Ji' S (<?) = J\\{q) - Jul 



(11) 



(12a) 



III. MEAN-FIELD CALCULATION 

We use a six-component vari- 

able S im = (ci m ,C2 m ,C3„ l ,si m ,S2„ l ,s 3m ). In a Fourier 
representation, the mean-field is 



j C 2 S (q) = J\\(q) + 



J- 



1:2 



Ji2 S +8J 2 c 3 s (<zr , (12b) 



# («) = Ml) + 3 J ™ + V J 1 C 2 + 8J 23 S (9) ■ (12C) 



3 



The eigenvectors are of the form 

v c 1 ' s (q) = (1,-1,0), 



(13a) 



respectively inducing a modulation in the columns, and 
favoring q c = 0. 



S (9)), (13b) 



(<?) = (j's" 5 - J\\{Q)> % 



J\\(q), 2Jg(q)). 



(13c) 



For each parameter set (J, G, J2), we numerically find 
which of the six eigenvalues is maximal and the corre- 
sponding q c . It turns out that the only two eigenvalues 
to be maximum are j| and jf . The critical temperature 
is thus the maximum of the following two temperatures: 



T c c - s = max(I# fl ) 
= i(cos(j, 

+||J±GU/l + 32cos 2 ig c . 



J 2 cos2 (Zc ) + |(J±G) 



(14) 



If T c c > T c s (resp. T c c < T c s ), only the cosine (resp. sine) 
components are involved in the transition. The bound- 
ary between the cosine and sine transitions is defined by 
T c c = T c s . Noticing the (a, a, b) structure of the u| and vj 
vectors, we conclude that columns 1 and 2 play similar 
roles, while column 3 has a distinct behavior. 

In the J-G plane, and for a particular value of J2, we 
identify four regions corresponding roughly to the four 
quadrants, two of them being shown in Fig. ^. The line 
G = 0, where the cosine and sine components are equiv- 
alent, is an obvious boundary. The curve 



\G\ = -g(J, J 3 ) J, 



(15) 



on which T c c = T', follows from a mechanism similar to 
the spin-flop mechanism of magnetism. For | J| and \G\ 
sufficiently high relative to |Ja|, q c vanishes, that is, the 
transverse couplings have "untwisted" the columns. The 
equation T c c (g c = 0) = T"(q c = 0) then gives a relation 
between J and G: 

(J + G) + \J + G\VS3 = {J - G) + \J - G\V33. (16) 



The solutions are G = and \G\ = — V33J. Thus, in 
this limit of q c = 0, g = V33 ~ 5.74, independently of J 
or J 2 . At q c — 0, the columns behave like a single vector 
flipping under the anisotropic effect of G. For smaller 
values of I J I and \G\, q c 7^ and g decreases, the same 
mechanism however remains. 

In the "cosine" regions, the helicity at the transition is 
a function q c = q c (J + G, J2), because J + G is the only 
combination of J and G appearing in jg. In the "sine" 
regions, we have q c = q c (J — G, J2) for the same reason. 
For J2 > — j, q c — 0, that is, each column stabilizes a 
"ferromagnetic" order p^Jl5[ . For J2 < — j, there is a 
region of the J-G plane where q c ^ 0, but q c = for \J\ 
and |G| high enough. The boundary is determined by the 
competition between \Ja\ and the transverse couplings, 




FIG. 2. Phase diagrams at the critical temperature for 
J-2 = -1.0, -0.8, -0.6 and -0.4. The full lines, all super- 
imposed with the precision used, are the boundaries between 
cos and sin phases and the dashed lines, the boundaries be- 
tween q c = and q c 7^ phases. For G < 0, the boundaries 
are a mirror image of the above with cos and sin phases being 
interverted. 



B. Thermal phase diagrams 

In order to investigate finite-temperature effects near 
the order-disorder transition, to better specify the nature 
of the phases, and to rule out the possibility of higher 
temperature first-order transitions, we set up a Landau 
theory from the microscopic model. We follow in essence 
the method proposed by Bak and von Boehm pjj . 



1. Free energy functional expansion 

For commodity, we divide the Landau free energy func- 
tional F into two parts: Ft and Fj. In reciprocal space, 
to fourth order, we find that 

f t =r^|^[c (U (gi)c A( (-gi) +s (U (gi)s /1 (-gi) 

+ \ 51 [ c P-(l^ c ^(l2)c tl (q 3 )c tl (~qi - q 2 - q 3 ) 
qi 92 93 

+2c f _ l (q 1 )c^{q 2 )s f _ l (q 3 )s f _ l (-q 1 - q 2 - q 3 ) 
+M?i) s A'( l ?2)s AI (<?3)s Al (-<7i - 12 - q 3 ) } (17) 

with T — 1/(3. C/j,(q) and s M (g) are the Fourier trans- 
forms of the mean values {c^ m ) and (s Mm ). This limited 
power expansion is numerically close to (less than 1% 



difference) the exact value up to v/(c^ m / 



^ 2 + (W 2 w 

0.5. A 6th-order development is 1% accurate up to 



y/ic^m) 2 + (s^m) 2 ~ 0.65. For every mean value un- 
der this limit of validity, a negligible number of ficti- 
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tious spins with modulus higher than 1 contributes to 
the statistics. We also find that 

Fj =~X]Z][ J ^(^)cm(9i)c I /(-9i) 



(18) 



These expressions represent the free energy functional 
for a group of three columns. The q summations have 
been truncated of all the umklapp terms. These umklapp 
terms would have pinned the modulation to commensu- 
rate values. By ignoring them, we allow incommensurate 
phases to occupy the entire parameter space, leaving a 
space of measure zero to commensurate phases. In real 
systems a devil's staircase [[lTj is expected instead of the 
continuous q profile. 



2. Order parameters 

We then assume that, near the transition, c^(q) = 
and Sfj,(q) = (V/x) except for q = ±q. In other words, 
we concentrate on the first harmonic to appear in the 
modulated phases. This is valid at the transition but it is 
not excluded that higher harmonics may appear at lower 
temperatures, as secondary order parameters. c M (g) and 
Sfi,(q) are the x and y components of three polarization 
vectors 

S M (g) = c M (g)x + s M (g)y. (19) 

These are complex quantities that may be expressed as 

c M = |c M |e^ and s tl = \s^\e i <. (20) 

This choice of variables allows for any elliptical polariza- 
tion and relative global phase for each column. To sim- 
plify notation, we replace |c M | by c M and |s M | by and, 
to avoid any ambiguity, we make no use of the complex 
and s M anymore. In real space, 



(S, k )=^S,(q)e^+S;(q)e- 



cos(qz k + ip°)x + s p cos(qz k 



Vu)y, (2i) 



with Zk = k + 5<5 M 3, so that ip^ s are the global phases of 
the different columns \l at the z = level. 

The function to minimize is then F = Ft + Fj, with 



F T =Tyh(cl + S l) + l\6(ci + 4) 



and 



+4[2 + cos2(^-^)]c^]} 



(22) 



ft V 



(23) 



with couplings as previously defined. F is a function of 
c M , Sjl,, ip^ s and q that, at first, is numerically minimized. 
The reader should note that the permutation 1 «-> 2 in 
the indices leaves F unchanged, which reflects the equiv- 
alence of columns 1 and 2. 



3. Helicity patterns 

We numerically observe simple relationships between 
the phases (p^ and tp s ^ . These, in turn, lead to a simplified 
expression for the free energy functional. For G = 0, the 
cosine and sine components are equivalent and s M = c M . 
By numerically minimizing F, we obtain 



- 1 - 2 ' 



(24) 



so that the modulation appearing in the columns is cir- 
cularly polarized, with a helicity given by the sign on the 
right-hand side. This sign (cr M = ±) may differ from one 
fi value to another. We denote (cri, <72, 03) the helicity 
configuration of the three sublattices of columns. For 
G = 0, the only allowed helicity configuration is (+ + +) 

(or equivalently ( )). Depending on the sign of J, 

the relative global phases of each component from one 
column to another are or ±7r. For J > 0, the columns 
adopt a "ferromagnetic-like" arrangement: 



C.S 



C.S 

v 3 » 



(25) 



and for J < 0, the triangular geometry imposes a colinear 
"antiphase" : 



C.S 

'Pi 



¥2 =V?3 



±7T. 



(26) 



Instead of the ordinary "120° state" for an evenly frus- 
trated system, our system concentrates the frustration 
in the 1-2 bond, which has a lower coordination number 
than 1-3 and 2-3 bonds. In contrast to the results ob- 
tained at T = in ||, where every column was forced to 
adopt the same amplitude and a deformed "120° state" 
was achieved, we obtained here a fully colinear antiphase. 

This fundamental difference arises from the freedom of 
the above model to adopt different amplitudes of modu- 
lation for the different columns. This was not allowed in 
S. Using the following definitions: 



C3 c 3 
Vc = — = — 
ci c 2 



S3 s 3 
Vs = — = — , 
si s 2 



(27) 



we have presented on Fig. for G = and for different 
values of J, r\ c and r\ s as functions of the temperature. 
Because G — 0, r\ c = r\ s = rj and the curves are identical 
for each J. This ratio has a well-defined value only below 
the critical temperature T c . It is to be noticed that for all 
cases presented, 77 is always larger than unity. For J < 0, 
the behavior is even larger than for J > 0. As a conse- 
quence, the displaced columns show a larger amplitude 
for the modulated phases. 
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These phase relationships divide the J-G plane in four 
quadrants, as seen in Fig. [|. This diagram possesses 
an athermal character since these relative phases are the 
only variables having an influence on the sign of each 
term of F, cos being always positive. However, it is 
important to stress that the phases (and the phase rela- 
tionships) have a physical meaning only when their cor- 
responding amplitudes are nonzero. From these partial 
results, we may rewrite F in a simpler form: 

F T = 2T{2c\ + 4 + 2s\ + s\ 

+ ±[3(2cf + 4 + 2s\ + si) + 2(2clsl + 44)]}, (32) 



FIG. 3. rj(T) curves for G = and different values of J. 
r/(T) is denned only for T below the critical temperature T c , 
depending on the J value. 

For G / 0, the rotational invariance is broken and we 
expect non-circularly polarized phases. Numerically, we 
still found tp^-ip^ = ±| and, if \ J\ > \G\, a (+ '+ +) 
configuration is realized. If | J| < |G|, we have a (+ H — ) 
configuration, as previously found in jll| (see Fig. |^). 
The inter-column relative phases are related to the sign 
of J + G for the cosines and of J — G for the sines. For 
J + G > 0, 

<Pt = <P°2 = V>3: (28) 
while for J + G < 0, 

tf=<fi=<pl± tt. (29) 

For J - G > 0, 

Vl = f2 = <Pl (30) 

while for J - G < 0, 

ifl = = ?l ± 7T. (31) 



G 



I 

+ -) 

±71) 



(+++) \ 


/ (+ + +) 


(±7C,±7t) ? 


\ (0, 0) * 



J 



(++-) 

(±71, 0) 



FIG. 4. Helicity patterns and phase differences in the J-G 
plane. The first number (0 or ±7r) is the difference ip%, — ipl, 
and the second is (fz — (fit . 



Fj = -(cosq + J 2 cos2q)(2cl + 4 + 2sl + s l) 

— 3( J + G)4 - 12\J+ G|(cos ig)cic 3 (33) 
-3( J - G)s\ - \2\J - G|(cos \q)s x s 3 . 

The absolute values of some couplings reflects the helicity 
choice made in order to minimize F. The positive signs 
of the coefficients of cfsf and 4 S 3 favorize competition 
between c and s variables. 



4- Phase boundaries at T c 

In this section, we investigate the first order-disorder 
transition to appear when lowering the temperature, that 
is, the highest temperature for which at least one order 
parameter is nonzero. This resulting critical temperature 
and the boundary between different phases should agree 
with the results found in Section A. This is presented 
in order to express the consistency of the Landau the- 
ory and its numerical treatment. For each parameter set 
(T, J, G, J2), we minimize F defined with the expressions 
d32| ) and (|33|); using (^) and fl23| ) would be equivalent, 
but numerically inefficient. The critical temperature is 
found by selecting T such that one of ci, C3, si or s 3 
is as small as possible, but nonzero. We have recovered 
the diagram of Fig. @: T c and q c being identical to those 
found in Section A. We also conclude that the transitions 
are always of second order, even with sixth-order terms 
in the expansion of Ft- 

For q c — 0, that is, when | J| or |G| are large enough, 
some analytical results as easily obtainable. The matrix 
{Hp}ij = {d 2 F/dxidxj}, where x is a vector constructed 
from the order parameters (x = (ci, C3, si, S3)), defines 
the local convexity of F. We diagonalize Hp to express 
this convexity along some eigendirections: hp, the four 
eigenvalues of Hp, are the convexity coefficients. For 
each parameter set (T, J, G, J2), F is minimal for a par- 
ticuliar x. For the algebric form of the present F, the sign 
of each convexity coefficient at x = indicates whether 
is a minimum or not: F is minimal at if every coefficient 
is positive, but one or more negative coefficients means 
that is no longer a minimum. A vanishing eigenvalue 
of Hp at x = corresponds to a second-order transition. 



G 



Hp has two distinct, doubly degenerate eigenvalues 
given by 

h c / =3[2T- (1 + J 2 ) - (J±G)] 

+ v/[2T - (1 + J 2 ) - 3( J ± G)] 2 + 144( J ± G) 2 . (34) 

ft,^ correspond to J + G and h s F to J — G. For a set 
(J, G, J2), T is the critical temperature T c when every 
hp is positive at x = 0, but at least one vanishes. The 
eigenvalues vanish at the following temperatures: 



rj-ic 



|(1 + J 2 



\(J±G) + || J ±G|V33. (35) 



We kept the solution that keeps T c positive for any J and 
G. These expressions coincide with ( |l4] ) if g c = 0. The 
cosine or sine components order depending on which of 
the two temperatures T£ or T c s is maximum. We have 
a phase boundary between the two types of order for 
T c c = T c s . This equality is equivalent to @. For q c ^ 0, 
the same analysis leads easily to (|l4|). 

For a small q c , we expand T^ ,s in powers of q 2 c . The 
coefficient is negative. The boundary between q c = and 
q c 7^ phases occurs when the g 2 coefficient vanishes. We 
have the relation 



22 



J±G\ = — (-J2-I), 



33 



(36) 



which defines the boundaries on Fig. B (the dashed lines). 



5. Thermal phase diagrams 

Below the critical temperature, and for J and G ex- 
actly lying on the |G| = —g(J,J2)J curve, both the c 
and s components order. For the same G and J 2 but 
for a higher J, the c components order first, followed at 
lower temperature by the s components, and conversely 
for a lower J. Consequently, the concomitant ordering 
of both the c and s components exists for a range of J. 
From the equivalence of the /i = 1 and = 2 columns, we 
have observed that c\ — c 2 and s\ — s 2 . Conversely, but 
because of their mutual hindrance, c\ < C3 and si < S3. 

We show in Fig. || the thermal phase diagrams for con- 
stant values of G and J 2 = — 1- The symmetry of F with 
respect to G <-> — G allows us to consider only G > 0. 
Each diagram has a high-temperature disordered phase 
(denoted d) and, at lower temperatures, ordered phases 
sin and cos. For q ^ 0, these phases are denoted Msin 
and Mcos. Just below the |G| = — g(J, J 2 )J curve is a 
mixed sm+cos phase that transforms continuously from 
sin to cos when approaching this region from large posi- 
tive and negative J. For q 7^ 0, the region sin+cos is an 
elliptical phase denoted E. 



a) 



1 5 



10 - 



b) 



10 




c) 



10 



G=l 




FIG. 5. Thermal phase diagrams for strong transverse cou- 
plings (see text). For G = 4 (a) there is no modulated phases. 
For G — 3 (b) the modulated phases are reentrant. The dot- 
ted lines indicate the limit of validity of the free energy limited 
power series expansion, below which at least one order param- 
eter is greater than 0.5. The insets schematically illustrate the 
orientations of the three molecules in a plaquette for sin and 
cos phases. 

For high values of |J| or |G|, q remains zero below 
T c , as shown in Fig. For small |G| and increasing | J|, 
steep boundaries between modulated and non-modulated 
phases are crossed (see Fig. ||c): q vanishes continuously 
from Mcos to cos and from Msin to sin. Inversely, for 
small |J| and increasing |G|, a reentrant boundary is 
crossed: a modulated elliptical phase reappears with in- 
creasing temperatures (see Fig. |]b). For high values of 
\J\, below the critical temperature, the disordered vari- 
ables (c in the sin phase and s in the cos phase) or- 
der at lower temperatures, forming sin+cos phases. For 
G = 1 (Fig. |c), we see the beginning of this phase for 
I J| >^ 2. The phase boundaries are parallel to the sin 
and cos boundaries with the disordered phase. These low 
temperature phases are not shown on Figs ^a and |5|b, but 
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they occur respectively for \J\ >~ 6 and 8. Because he- 
licity configurations have meaning only when c M , s M and q 
are nonzero, that is, for E phases, the diagrams of Fig. || 
do not show the richness of Fig. ||. For the choice of 
parameters, the central E phases in Fig. [| are (+ H — ) 
phases. For smaller values of \G\, it is expected that a 
phase boundary between E phases having (+ H — ) and 
(+ + +) configurations would be observable. 

Indeed, the quasi one-dimensional regime, where 
\Ji\, ^> \J\,\G\, is interesting because it is more re- 
alistic for HHTT, with E phases having different helicity 
patterns. We show in Fig. ^ thermal phase diagrams for 
G from 0.1 to 0.5. Both the E{+ + +) and E{+ + -) 
phases exist and are sometimes adjacent. On the J inter- 
val presented in Fig. ||, all the ordered phases are mod- 
ulated. The temperature "depth" of the Msin and Mcos 
phases is approximately proportional to \G\ and indepen- 
dent of J. 



3 
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FIG. 6. Thermal phase diagrams in the quasi-onedimen- 
sional regime (see text). For G = 0.5 (a) the E(+ H — ) and 
E(+ + +) are separated. For G = 0.3 (b) and G = 0.1 (c), 
dashed lines denote first-order helicity reversal transitions. 

Fig. [7] illustrates schematically the positions and orien- 
tations of the molecules for different phases encountered 
in the diagrams of Fig. 0, in the case where G > 0. Some 
are also present in Fig.^j. The length of the tails repre- 
sent the amplitude of ordering. The Msin phase (Fig. 0a) 
shows modulated order with uniform orientations along 



the columns and where column 3 is in phase opposition. 
For G > 0, the Mcos phase (Fig. 0b) is a modulated 
ferromagnetic-like state. The maximum lengths of the 
tails shows that column 3 is more ordered than columns 
1 and 2. The corresponding sin and cos phases are il- 
lustrated in Fig. ||a. Fig. (t]c to d represent the various 
elliptical phases E encountered from left to right in Fig. @. 




FIG. 7. Schematic illustration of the orientations of 
molecules and amplitudes of mean values with G > for 
a) the Msin phase with column 3 in phase opposition, b) 
the Mcos phase with column 3 in phase conjunction, c) the 
E(+ + +)(tt,ty) phase, d) the E(+-\ — )(0,7r) phase and e) the 
E(+ + +) (0,0) phase. 



IV. MONTE CARLO SIMULATIONS 
A. Method and algorithm 

We use the Metropolis algorithm (see |Q on general 
Monte Carlo methods) to simulate a three-dimensional 
lattice of XY variables. Our goal is to obtain the essential 
features of the thermal phase diagram of the system (0), 
restricting ourselves to the quasi one-dimensional regime. 

The above method has been implemented in conjunc- 
tion with the "spiraling" algorithm [^,^0|, which simu- 
lates incommensurate helical phases by relaxing the con- 
straint associated with the finite size of the lattice. Out- 
side the critical regime, the results are in principle in- 
dependent of the lattice size. This method is applied 
on each column of N XY variables {cf>i, ...,</> at}. We 
simulate the neighbors of (f>\ and 4>i with ■phantom sites 
(denoted by a prime) related to sites of the opposite ex- 
tremity of the column: 
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N—l 



m-i 



NA, 



J N 



NA. 



(37) 



For the sites <f>N—i an( l 4>n, we use the phantom neighbors 



■NA, 



o 2 = q> 2 



NA. 



(38) 



The reader will note that, for A = 0, simple periodic 
boundary conditions are recovered. A is an effective field 
representing an additional indefinite length of a lattice 
modulated with a constant pitch. For the system to se- 
lect its own boundary conditions, we consider A as a 
thermodynamic variable. Each Monte Carlo step of the 
spiraling algorithm consists in N ordinary trial flips of 
the <j) variables and a single trial flip of A. We modify A 
by a small random angle 5 A. In order for the new A to 
be compatible with the above equations, a twist must be 
imposed on the lattice: 



4>, <- (f>i + (i - 1)5A. 



(39) 



The new state has a new total energy, and is then tested 
for acceptance with the Metropolis algorithm. We may 
compare the spiraling algorithm to a high-order mean- 
field approximation, where the exact clusters have the 
size of the finite lattice used in the simulation. Close to 
a transition, when the correlation length is large relative 
to the cluster size, this approximation loses its validity 

We simulate a 6 x 6 triangular lattice of columns, that 
is, 12 groups of three columns. We use periodic boundary 
conditions in the plane and the spiraling algorithm in the 
columns' direction. Each column has its own A variable. 

For two neighboring columns with non-planar inter- 
actions ( J' and G"), the spiraling algorithm introduces 
some small energy discrepancies. This problem arises 
from the non-equivalence of some of the couplings at the 
edges of the lattice. Due to the relatively small energy 
involved in this boundary effect (especially for a quasi 
one-dimensional system), it was neglected. 

We identify the different phases with the nonzero val- 
ues of the Fourier coefficients of the x and y components 
defined by: 



m) 



12N X X/ 



COS (pijk COSqZfc, 



(40a) 



12N ^ ^ YlJ 

(j,i)efl„ k 



cos 4>ijk sin qzk, (40b) 



Using these coefficients, we construct the order parame- 
ters: 



(41a) 



(41b) 



The (c^(q)} and (s^q)) profiles show a maximum at 
q = qo, and may show secondary peaks. Indeed, it has 
been shown |2l[] that a third harmonic at q% = Zq® should 
appear for linearly polarized columns. In our simulations, 
this third harmonic should be an indication of such a lin- 
ear polarization. We also construct the relative angular 
phases of the columns defined by: 



cos 



(&-&)(?) 



<(gK(g) + feffi)*£(g) 



(42a) 



Bm(0J-^)(«) = — 



s^c^q) 



(42b) 



cos! 



<t>im 



<+M<{ql + K+i(q)b c M 



c^+iiqy^q) 



(42c) 



<V+i 



<V+i 



<P S M) 



^+l(9K(9) -a£+l(9>£(9) 



c M +i(9)c M (g) 



(42d) 



<+i(q)am+bl +1 (q)b s Jq) 



(42e) 



sm 



>+i 



KM 



^+1(9)^(9) -0^+1(9)^(9) 



1(9)5^(9) 



(42f) 



b s — 0„)(9o)) to the mean-field expres- 
and so on for each mean value. 



127V 



We compare (cos( 
sion cos(<£* -(/>£). 

The amplitudes (41) converge relatively fast, but the 
angular phases ([i"2] ) converge much more slowly: they are 
ratios of fluctuating quantities. For some simulations, 

val- 



sin <f)jjk cos qZk , (40c) we are only interested in the (c^(q)) and (s^(q)} 



12N 



sm (pij k sm qz k , 



(40d) 



with Zk = k + ^6^3. i? M represents the sublattice of /1 
columns and N is the number of sites in each column. 



ues, and we have simulated only 2000 MCS/S, including 
1000 MCS/S for thermalization. N, the number of sites 
per column, was taken to be 40. However, many more 
steps was needed to obtain the relative angular phases 
(25000 MCS/S, including 5000 MCS/S to thermalize). 
In these cases, iV was taken to be 12. To improve numer- 
ical efficiency, we discretize the (j) an d A variables into 
256 values from to 2n. 
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All simulations are done with G = 0.1 and J2 = — 1, a 
quasi one-dimensional limit. We perform the simulations 
by gradually decreasing the temperature for each value 
of J. This allows a greater numerical stability for low- 
temperature phases and retains the helicity sign ((+ + = F) 
or ( ±)) for all temperatures. 

Preliminary simulations on a single plaquette of three 
columns were done to compare the spiraling algorithm 
with periodic boundaries conditions in some reasonable 
computing time. For N — 40, simulations using periodic 
boundaries along columns give very similar results com- 
pared to simulations with the spiraling algorithm, except 
for smaller mean values (because of higher fluctuations) 
and slightly displaced peaks. However, the difference is 
pronounced for N = 12, where the spiraling algorithm 
broadens the principal peak, while periodic boundaries 
conditions destroy the whole spectrum. 




FIG. 8. a) (ci(g))(T) and b) (si(q))(T) profiles for 
J = —0.15, G = 0.1 and J2 = — 1. The arrows indicate the 
first and third harmonics. The qo peaks appear at T w 0.4 
for ci and T ~ 0.9 for si. 



B. Thermal phase diagrams 

Fig. H shows a typical result for the order parameters, 
from which the thermal phase diagram is constructed: 
the profiles (a(q))(T) and (si(g))(T). Notice that both 
go and (73 peaks appear at specific temperatures ((73 is 
folded in the [0, 7r] interval). The (c M (g)) and (s M (g)) 
profiles show that, for J = —0.15, a decreasing tempera- 
ture drives the system from d phase, to Msin phase and 
finally to E phase. The critical temperatures are arbi- 
trarily taken to be the points at which the amplitude of 
the peak is half its maximum value. 

The thermal phase diagram is constructed by repeat- 
ing the simulations for many J values. Fig. ^ shows the 
diagram for G = 0.1 and J 2 = — 1. We always obtain 
(cos(^» -^)(« )) » and (sin(^-^)(g )) « ±1, with 
signs corresponding to the mean-field helicity configura- 
tions: (+ + -) or ( h) for I J\ < 0.1 and (+ + +) or 

( ) for I J I > 0.1. Some longer simulations were done 

to establish clearly the phase differences, at J = —0.15, 
—0.05, 0.05 and 0.15. Given some high-enough (c M (<?o)) 
and (s/i(?o)) amplitudes, every relative phase is compat- 
ible with the mean-field calculations and the results just 
cited. Moreover, many short simulations where done on a 
single plaquette of three columns (N = 40) and, although 
no phase transition is observed (the system is onedimen- 
sional), the phase differences arc in excellent agreement 
with Fig. |. 

Contrarily to the mean-field results (Fig. |6|c), the E 
phases with different helicity configurations were not 
found to be adjacent: a linearly polarized phase opens 
up between the E{+ + +) and E(+ H — ) phases. On the 
diagram, for a given value of J, a dot at T — means that 
no phase transition was clearly identified when lowering 
the temperature. 
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FIG. 9. Thermal phase diagram for G = 0.1 and J2 = — 1 
from Monte Carlo simulations. The continuous line is a guide 
to the eye (see Fig. for the mean- field diagram). 

V. DISCUSSION AND CONCLUSION 

X-ray measurements on HHTT have lead to the obser- 
vation of a sequence of phases as temperature is lowered. 
Two neighboring phases {Dhd and H) at intermediate 
temperatures have a triangular structure of columns of 
molecules. Entering the H phase from the Dhd phase in- 
volves a mutual and concomitant ordering of both the po- 
sition and orientation of the molecules along the columns. 
However, it is immediately realized that the positions of 
the molecules are rapidly frozen as compared to their ori- 
entations, being submitted to more stringent intermolec- 
ular forces. This is the main justification for the model 
studied in this paper, limited to orientational degrees of 
freedom. The obtained results and in particular the gen- 
eral trend of the thermal phase diagram should be useful 
to understand the behavior of HHTT inside the H phase. 

The first result of interest is the ordering of linearly 
polarized phases at T c with a finite wavenumber for the 
amplitude modulation in the columnar direction, whose 
value is decreasing with increasing values of the inter- 
molecular couplings and eventually vanishes. A similar 
disappearance of the amplitude modulation has been pre- 
dicted Q at T = where it then shows up as an unwind- 
ing of the helical pitch. At T c , for large \ J\ and |G|, a 
boundary at a constant slope of magnitude ^/33 is pre- 
dicted between two linearly polarized phases. It is to be 
noticed that this behavior may be seen as a reminiscent 
effect of the two boundary structures predicted at T = 
between linearly polarized phases. At T = 0, the con- 
stant slopes of the boundaries are respectively —3 and 
—5. As seen on Fig. |^a, this T = behavior is perfectly 
predictible at large G. 

The frustration between ordered columns of molecules 
with an octupolar moment on a triangular lattice is very 
high. Part of this frustration is relaxed by freezing the 
positions of the molecules and displacing one of every 
three columns by half a lattice spacing in the columnar 
direction. However, for negative value of the interaction 
parameter J, even under these conditions, substantial 
orientational frustration remains. At T = 0, with only 
constant amplitude phases, a non-colinear distorted 120° 
phase was obtained with global relative angular phases 
between the columns which depended on the pitch of 
the helical modulation of the columns. In the present 
case, allowing for different amplitudes of the modulation 
for the displaced columns compared to the undisplaced 
columns, the resulting configuration is colinear and the 
relative angular phases is independent of the pitch. This 
angular configuration, typical of unfrustrated systems, is 
only achieved through a larger amplitude of modulation 
for the displaced columns as compared to the undisplaced 
ones, as shown if Fig. 0. 



For moderate values of G as compared to the intra- 
columnar couplings (Fig. ||), when there exists elliptical 
phases E, both the helicity patterns (+ + +) and (+ + — ) 
are predicted. Also it is to be noted that the temperature 
range over which the linearly polarized and modulated 
phases exist, narrows down with decreasing value of G. 
Under these conditions, it is expected that we rapidly 
enter the elliptical phases on lowering the temperature 
below T c . In this temperature range, where our model, 
which freezes the positional degrees of freedom rapidly 
under T c , is certainly valid, and for \G\ > |J|, the ellipti- 
cal phase (+ H — ) is predicted. Recall that the (+ H — ) 
phase is the one observed Q] for HHTT in the H phase. 
We then conclude that HHTT is a quasi one-dimensional 
system where the molecular orientations are dominated 
by the non-rotationally invariant interaction between oc- 
tupolar moments located on a distorted triangular lat- 
tice. 

The finite temperature Monte Carlo simulations on fi- 
nite size systems has confirmed remarkably the mean field 
results. The one important difference is that the (-M — ) 
and (+ + +) configurations are not seen to be adjacent 
in Monte Carlo simulations where they were predicted 
to have a common boundary in the mean field approx- 
imation. The hard-spin constraint |S Mm | = 1 is auto- 
matically satisfied in the Monte Carlo simulations. At 
low temperature and in the mixed regime |J| ~ |G|, 
the system is indecisive between configurations where all 
the columns have identical helicities and opposite hclici- 
ties between the displaced columns and the undisplaced 
columns. Based on the hard-spin contributions which are 
optimal for non-helical phases, the system seems to be 
driven to zero temperature while retaining either a mod- 
ulated cos or sin phase according to the sign of J. This 
result probably manifest itself into the rare occurence of 
helical phases experimentally. 

Regarding the Dhd — * H transition in HHTT, this 
work raises the question of the detailed nature of the 
observed H phase near the transition: is there a linearly 
polarized Ll H phase"? The presence of a linearly po- 
larized modulated phase would be an indication of the 
chiral octupolar nature of the molecule (the G coupling). 
Such phases could be observed as a third harmonic in the 
modulation of the columns in X-ray measurements. The 
most important result of this paper is the higher ampli- 
tude predicted for the displaced columns as compared to 
the undisplaced columns. This predicted behavior would 
affect the X-rays results and in particular the relative 
amplitude and thermal broadening of the Bragg's peaks. 
The J and G values may be varied by using discoid com- 
pounds other than HHTT or by the application of surface 
pression. 



11 



ACKNOWLEDGMENTS 



The authors are grateful to M.L. Plumer for many dis- 
cussions. G.L. would like to thank W.M. Saslow for help- 
ful suggestions related to the simulations. This work is 
supported by NSERC (Canada) and by FCAR (Quebec). 



[1] S. Chandrasckhar, Liquids crystals (2nd edition), Cam- 
bridge University Press (1992). 

[2] P.-G. de Gennes and J. Prost, The physics of liquid crys- 
tals (2nd edition), Oxford University Press (1993). 

[3] E. F. Gramsbergen, H. J. Hoving, W. H. de Jeu, K. Prae- 
fke and B. Kohne, Liq. Cryst. 1, 397 (1986). 

[4] E. Fontes, P. A. Heiney and W. H. de Jeu, Phys. Rev. 
Lett. 61, 1202 (1988). 

[5] E. Fontes, Ph.D. dissertation, University of Pennsylvania 
(1989). 

[6] S. Idziak, Ph.D. dissertation, University of Pennsylvania 
(1992). 

[7] A. Caille and M. Hebert, Phys. Rev. E 54, R4544 (1996). 
[8] M. Hebert and A. Caille, Phys. Rev. E 53, 1714 (1996). 
[9] M. Hebert and A. Caille, Phys. Rev. E 51, R1651 (1995). 
[10] M. Hebert and M. L. Plumer, Phys. Rev. E 54, 550 
(1996). 

[11] M. L. Plumer, A. Caille and O. Heinonen, Phys. Rev. B 

47, 8479 (1993). 
[12] M. Pesquer, M. Cotrait, P. Marsau, V. Volpilhac, J. Phys. 

France 41, 1039 (1980). 
[13] M. Cotrait, P. Marsau, M. Pesquer, V. Volpilhac, J. Phys. 

France 43, 355 (1982). 
[14] W. Selke, Spatially Modulated Structures in Systems 

with Competing Interactions, in Phase Transitions and 

Critical Phenomena, vol. 15, Academic Press (1992). 
[15] S. Redner and H. E. Stanley, Phys. Rev. B 16, 4901 

(1977). 

[16] H.-Y. Choi, A. B. Harris and E. J. Mele, Phys. Rev. B 

40, 3766 (1989). 
[17] P. Bak and J. von Boehm, Phys. Rev. B 21, 5207 (1980). 
[18] K. Binder (editor), Monte Carlo Methods in Statistical 

Physics, Topics in Current Physics 7 (1979). 
[19] W. M. Saslow, M. Gabay and W.-M. Zhang, Phys. Rev. 

Lett. 68, 3627 (1992). 
[20] M. Collins and W. M. Saslow, Phys. Rev. B 53, 8533 

(1996). 

[21] M. B. Walker, Phys. Rev. B 22, 1338 (1980). 



12 



